Mesoscopic photon heat transistor 
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We show that the heat transport between two bodies, mediated by electromagnetic fluctuations, 
can be controlled with an intermediate quantum circuit - leading to the device concept Mesoscopic 
Photon Heat Transistor (MPHT). Our theoretical analysis is based on a novel Meir-Wingreen- 
Landauer type of conductance formula, which gives the photonic heat current through an arbitrary 
circuit element coupled to two dissipative reservoirs at finite temperatures. As an illustration we 
present an exact solution for the case when the intermediate circuit can be described as an elec- 
tromagnetic resonator. We discuss in detail how the MPHT can be implemented experimentally in 
terms of a flux-controlled SQUID circuit. 



Problems involving interactions between quantum par- 
ticles and electromagnetic fields have a long and rich 
history. During the last decades developments in meso- 
scopic physics have provided unprecedented possibilities 
to engineer the electron-photon interactions. A major ad- 
vantage in mesoscopic systems is their versatility which 
allows a high-level of control in their design and oper- 
ation. Recent advances in this field include impressive 
quantum state manipulations involving single photons in 
circuit cavity QED experiments [ij and a demonstration 
of single-channel photon heat transport [2| . 

At low temperatures, when phonon modes become ef- 
fectively frozen, photonic heat conduction becomes the 
dominant channel for thermal transport @> I n this 
Letter we study photonic heat transport in a structure 
consisting of two reservoirs coupled via an intermediate 
electric circuit. The reservoirs are assumed to behave as 
linear dissipative circuit elements and are thereby fully 
characterized by their response functions and tempera- 
tures. Furthermore, at low temperatures the wavelengths 
of relevant field fluctuations are much longer than a typ- 
ical system size so the reservoirs can be effectively con- 
sidered as lumped elements. Within these assumptions 
we apply the Caldeira-Leggett procedure and model the 
reservoirs as continuous distributions of harmonic oscilla- 
tors. By applying nonequilibrium Green's function meth- 
ods we derive a formally exact Meir-Wingreen-Landauer- 
type formula 0] for the heat current through the struc- 
ture. Our formula involves the noise power of the inter- 
mediate circuit, in the presence of the coupling to the 
leads, and serves as a general starting point for solving 
the heat transport problem. We find an exact solution 
for the heat current flowing through an electromagnetic 
resonator circuit and show explicitly that, in analogy 
with the transistor effect in charge transport problems, 
the heat flow through the structure can be modulated 
by applying an external control to the middle circuit. 
We suggest that an experimental demonstration of the 



heat-transistor action can be achieved by using a Super- 
conducting QUantum Interference Device (SQUID) cir- 
cuit as the tunable resonator. Thus, the magnetic flux 
controlled heat current is a photonic analogue to the 
gate voltage controlled electronic heat current recently 
demonstrated in Ref. [5|. 

Now we turn to the technical derivation of the general 
formula for the heat current in the system consisting of 
a left and a right reservoir, and an arbitrary quantum 
circuit between them (Fig. [I}. We treat the problem 
by employing a nonequilibrium Green's function method 
analogous to those used earlier in electron and heat trans- 
port problems [1, 0, B 0]- The reservoirs are de- 
scribed by quadratic boson fields and, according to the 
Caldeira-Leggett prescription, can be thought of as ar- 
bitrary linear electric circuits by choosing specific distri- 
butions of frequencies u)j and couplings gj (introduced 
below) 13, U | • The total Hamiltonian is assumed to be 
of the form H = Hl + Hr + Hm + He , where 



Htjr — 



1/2), 



and the inductive coupling term is 
H c = MI (i L + i 



ill) , 



(1) 
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which involves the current operators for the central re- 
gion I and for the reservoirs Il/r = J2jeL/R9j(^j + ®j)> 
respectively. The specific form of the Hamiltonian and 
the current operator of the middle region do not need to 
be specified at this point (below we shall treat a specific 
example). The mutual inductances between the middle 
circuit and the leads are assumed to be equal, though 
this is not necessary in the following derivation. A ca- 
pacitive coupling between the reservoirs and the middle 
system can be treated in close analogy with the inductive 
coupling studied here [121 ]. 



The reservoir Hamiltonians Ht,/r do not commute 
with the total Hamiltonian H, thus giving a rise to an 
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FIG. 1: Temperature gradient between the left and right 
reservoirs induces photonic heat current through an arbitrary 
quantum circuit coupled to them. Heat flows from left to 
right when Tl > Tr. The reservoirs are assumed to behave 
as linear dissipative elements and they couple only through 
the middle circuit. In principle a direct coupling between the 
reservoirs always exists but in practice it can be made negli- 
gible by an appropriate sample design. 



energy flow in the structure. This energy flow is charac- 
terized by a heat current Jl / R defined as 



J L / R (t) = (Hl/r) = iM Y, [9j^{aj(t)I(t)) - h.c 
jeL/R 

= -2MRe Y 9jUjGf(t,t), 



(3) 



jeL/R 



where G">(t,t') = —i{a,j(t')I(t)) . The transport problem 
is reduced to finding the lesser Green's function G< (t, t') 
which can be derived by the equation-of-motion tech- 
nique @, Il3j |. Following the standard prescription Q, 
we first consider the equilibrium zero-temperature time- 
ordered correlation functions, which obey 



{id v - u 3 ){T[a 3 {t')I{t)]) 



M 9j 



(T[I(t')I(t)}), (4) 



or, after a formal integration, 
Mg 3 



(T[a 3 (t')I(t)}) 



dti{T[I(t)I{h)])D j (ti--e), 



(5) 

where Dj{t\ — t') is the free reservoir Green's function. In 
nonequilibrium, this equation holds on the Keldysh con- 
tour, and using the analytical continuation rules known 
as Langreth's theorem [fj], we obtain 



Gf(t,t') = 



M 9j 



dt x 



(i^iit^YDf^-t'y 



(7(f)/(t 1 ))< J D«(i 1 -0 



(0) 



where the superscripts r, a and < stand for "retarded", 
"advanced" and "lesser", respectively. Explicitly, the 
current correlation functions are {I(t)I(t')) r = —i0(t — 
t')([I(t)J(t')]) and (i(t)I(t / ))< = -i</(f )/(<)). In a 
steady state Gf(t,t') = Gf(t — t'), and it is convenient 
to introduce the Fourier transform: 



Gf{u>) 



M 9j 



{II) r (w)Df(w) + (jl)<(w)m(w) , (7) 



where D°-(lo) — l/(ui—uij—ii]), D<(lo) = —i2irn(ujj)8(uj- 
ujj), and n(uj) is the Bose function. We thus obtain 



2nh 
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-Im(7 1) r ' {uj)2i:ri{uj j)S(uj — ujj) 



+lm(ii) < (uj)7rS(uj - Uj) 



(8) 



To make a connection to quantities with a clear physical 
interpretation, it is useful to express Eq. ([8]) in terms of 
the noise power. The noise power for an observable A 
is defined as S A (iu) = j '^dt e iu - '(*-*') (A{t)A{t')) , so the 
current noise in the left lead is 



+(n(uij) + 1)2tt(5(cj - uij)} 
In terms of the noise power Eq. ([8]) becomes 
duju) 



(9) 



j L = 2Ar 
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-lm(Iinco)S lL (-io) 
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+1m(II)<(u)-(S iL (u)-Si L (-w)) 



(10) 



Also the /-correlation functions can be written in terms 
of noise power: — lm(I I) r (ui) = |(S/(w) — Si(-uj)) and 
lm(II) < (u>) = —Si(-uj), which yields 



Jl = M 



duuj 



1 ^ ^| [Si(u)S iL (-u) - 5x(-w)5 i£ (a;)] . 



(11) 

A similar expression holds for J R ; in steady-state situ- 
ations J = Jl = —Jr.- Formulas (fTTJ)) and (fTTj) con- 
tain the free reservoir noise functions which can be ob- 
tained straightforwardly from the admittances Y l / r (lj) 
and temperatures Ttjr by applying the Fluctuation- 
Dissipation theorem [lfj: 

S iL/R (u) = Re[Y L/R (u>)]hu> [coth(/3 L/J? ^/2) + l] . 

(12) 

With the help of Eq. (fl"2"|) , the heat current (fTl) can be 
written as 

Jl = I — {2 [S r (u>) - Sj(-w)] Re[Y L (u,)]n L (u;) 



2tt 
-Si(-w)2Re\Y L (u)]} 



(13) 



If the lead admittances share the same frequency depen- 
dence, Yl(oj) = cY R (ui) with some constant c, the sta- 
tionary heat current can be expressed as J = Jl/(c + 
1) — cJr/ (c + 1) and cast into the Landauer form: 



J =M 2 



diouj 2 
2tt 



[S x {u) - Si(-w)] 



2Rc[Y L (m)}Re[Y R (u)} 
Be[Y L {u))+Re\y R (uj)] 



[ti l {uj) - n R (w)]. (14) 
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Equation ([14]) is the main formal result of this Let- 
ter. We emphasize the strong analogue between this re- 
sult and the Meir-Wingreen conductance formula: here 
the "bosonic thermal window" (jil — n R ) replaces the 
"fermionic voltage window" (n£ — n^), the real part of 
the admittance plays the role of the line-width function 
r ' L/R-, an d the noise power Si(u>) — Si(—u>) of the central 
region, which needs to be evaluated separately, replaces 
the central region spectral function. Since Si contains in- 
formation of the internal dynamics of the middle region 
in the presence of coupling to the reservoirs, its evalu- 
ation may be difficult indeed, and only in special cases 
analytic progress can be expected. In the following we 
calculate Si(uj) for an electromagnetic resonator and il- 
lustrate how the reservoir admittances come into play 
through the self-energy of the middle circuit. 

Suppose now that the mediating quantum circuit is an 
electromagnetic resonator with inductance L and capac- 
itance C (see Fig. [2] a). The central region Hamiltonian 
then takes the form Hm = fvu>o{lrb + |), and the current 
operator is / = Io(b+w), where Iq = y // hujo/2L and uq = 
1/y/LC; b, w are bosonic creation and annihilation oper- 
ators, [b,w] = 1. To evaluate p^|) one needs to find the 
retarded function (I(t)i(t')) r = -i8(t - t')([I(t)J(t')}) 
which can be expressed as sum of four different retarded 
functions 



The reservoir admittances are given by the Kubo formula 



(//») = 1% 



(btf) r {L0) + {Pb) r {uj) 



(15) 



Note the presence of the anomalous functions (b(t)b(t')) r 
and (b'{t)b'(t')) r ; they are important and come into 
play because the interaction term also includes the non- 
rotating- wave terms cub and Ojb^ . By using an equation- 
of-motion technique [y, [Hj we find a closed set of equa- 
tions 



(w - oj + i V )(btf) r = i + s r H (ffly + (bWy 
(LU + UJ0 + M ?)(& t & t ) r = -s r M [(SS t ) r + (tftfy 

(cj + uj + ir)){Pb) r = -1 - E r H \{bb) r + (b r b) r 
(uj-LO + in)(bb) r = E"» [(bb) r + (Vb) r ] , (16) 

where the retarded self-energy is given by 



(MIo) 2 



h 2 ~ J _ ' JJ \uj — ujj + if] to + ujj + ir\ 



i{iLih) r {^) 

— V 

huj J \u — ujj + in u> + uij + in 



(18) 



so the self-energy is related to the reservoir admittances 
through relation 

h 

The algebraic system (fTB| is readily solved yielding 

—ifko 



<//») = 



ujF(lo + in) + M 2 u 2 (Y L (u) + Y R (u)) ' 



(20) 

where we have introduced shorthand F(u>) = 
(ih/lQ)(oj 2 — Uq)/(2uiq). By extracting the imag- 
inary part of Eq. (|20|) and recalling the relation 
-2Im(//) r (w) = Si(u) - Si(-w) we find 

_ f°° duj Ahw 5 M i Re[Y L {oj)}Re[Y R {uj)] 

~ Jo 2^ \loF{uj) - M 2 u 2 {Y l {uj) + Y R {io))\ 2 * 

x [n L (uj) - n R (w)] . (21) 



The result (|2T|) is valid for arbitrary reservoirs admit- 
tances, as long as the intermediate circuit can be de- 
scribed as an oscillator. It can be shown that the upper 
limit of heat current (|2ip is given by the universal single- 
channel heat conductance quantum Gq — irk%T/3h, 



(17) 



0- 0, @, Ei , IE 3| , an exact parallel with the electrical 
single-channel conductance Go = 2e 2 /h, which follows 
from the Meir-Wingreen formula for a noninteracting res- 
onant level model at resonance [f|. 

In the following we assume that the reservoirs are iden- 
tical and consist of a resistor, a capacitor and an induc- 
tor in series. Thus the admittances are can be written 
as Y l (uj) = Y R {u) = 1+iQ{ ^ 1 r _ Ur/ujV where R, Q and 
lo r are an effective resistance, a quality factor and a res- 
onance frequency of the reservoir. The net heat flow is 
determined by five dimcnsionlcss parameters M 2 luq / LR, 
u>o/u) R , Q, hjjQ / koTL and hjjQ/k^T R where Tl and T R 
are temperatures of the reservoirs. By externally con- 
trolling the parameters of the middle circuit it is possible 
to tune the heat flow through the structure. A prac- 
tical way to realize this is to employ a dc-SQUID as a 
middle circuit and apply a magnetic flux <I> through it 
(see Fig. [U b)). At low temperatures the system can 
be modeled by an LC-oscillator with a tunable induc- 
tance L($) sa i/|cos(7r$/$o)| where $ = h/2e is the 
flux quantum and L is determined by the Ambegaokar- 
Baratoff critical current [ljj]. The LC-oscillator descrip- 
tion of the SQUID circuit is expected to be accurate 
within realistic parameter values at the experimentally 
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FIG. 2: An electromagnetic resonator as the intermediate cir- 
cuit (a)). In b) the mediate resonator circuit with a tunable 
inductance is realized by a dc-SQUID with an external mag- 
netic flux bias $. The external flux can be used to modulate 
the heat current through the structure and plays an analogous 
role to the gate voltage in electronic transistor. 



verified crossover temperature T cr ^100 mK below which 
the photonic thermal conductance should dominate Q. 
A numerical evaluation of Eq. pTj) is presented in Fig. [3J 
which shows the heat flow as a function of the external 
bias flux <£>. The maximum flow is obtained at at in- 
teger values of &/&o, whereas at half- integer values the 
reservoirs are thermally decoupled. At low quality fac- 
tors the reservoirs are more efficiently matched and the 
system has a better thermal coupling. By increasing the 
coupling parameter M 2 ujq/ LR the maximum value of the 
heat flow could be enhanced closer to the single-channel 
maximum value. 

In summary, we have studied photon heat transport in 
nanoelectronic circuits based on a novel Meir-Wingrccn- 
Landauer formula in a two-terminal geometry. This for- 
mula expresses the heat current in terms of the admit- 
tances of the heat reservoirs, and the noise power of the 
intermediated mesoscopic circuit. The formula can serve 
as a starting point for an analysis of photon heat trans- 
port in a wide range of applications. As an example, we 
present an exact solution to the transport problem in the 
case of an electromagnetic resonator playing the role of 
the mediating circuit. We propose a new device concept, 
a mesoscopic photon heat transistor, where the heat cur- 
rent through the structure can be strongly modulated by 
the external magnetic flux through a dc-SQUID loop. 

We would like thank Antti Niskanen, Jukka Pekola and 
Tero Heikkila for valuable discussions. 
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FIG. 3: Heat flow through the SQUID structure as a function 
of applied flux $, corresponding to parameters M 2 u>o/LR = 
1, ujq/ujr — 1 and Tr = Tl/2 = huo. The different curves cor- 
respond to different reservoir Q values, Q — (blue), Q = 0.1 
(red), Q = 0.5 (green), Q = 2 (black) and Q = 10 (cyan). The 
heat flux is normalized with respect to the universal single- 
channel maximum value J max = ^{Tl - T%). In the inset 
the heat flux is plotted as a function of the temperature of the 
right reservoir corresponding to the parameters uiq/ujr = 1, 
T L = 2hu , Q = 0.1 and $ = mod $ . The different 
curves correspond to M 2 u>o/LR = 0.5 (blue), M 2 ujq/ LR — 1 
(red), M 2 to /LR = 3 (yellow), M 2 lu /LR = 5 ( green) and 
M 2 lu /LR = 10 (black). 



